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The sensitivities to surface friction and the Coriolis parameter in tropical cyclogenesis 
are studied using an axisymmetric version of the Goddard cloud ensemble model. The 
Coriolis parameter is related to earth's rotation rate and the latitude. Our experiments 
demonstrate that tropical cyclogenesis can still occur without surface friction. However, 
the resulting tropical cyclone has very unrealistic structure. Surface friction plays an 
important role of giving the tropical cyclones their observed smaller size and diminished 
intensity. Sensitivity of the cyclogenesis process to surface friction, in terms of kinetic 
energy growth, has different signs in different phases of the tropical cyclone. Only in 
the brief initial onset phase can increase surface friction yields faster kinetic energy 
growth. Contrary to the previous theory of Ekman pumping efficiency, which implies a 
preference for the highest Coriolis parameter in the growth rate if all other parameters 
are unchanged, our experiments show no such preference. 
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Abstract 


The sensitivities to surface friction and the Coriolis parameter in tropical 
cyclogenesis are studied using an axisymmetric version of the Goddard cloud ensemble 
model. Our experiments demonstrate that tropical cyclogenesis can still occur without 
surface friction. However, the resulting tropical cyclone has very unrealistic structure. 
Surface friction plays an important role of giving the tropical cyclones their observed 
smaller size and diminished intensity. Sensitivity of the cyclogenesis process to surface 
friction, in terms of kinetic energy growth, has different signs in different phases of the 
tropical cyclone. Contrary to the notion of Ekman pumping efficiency, which implies a 
preference for the highest Coriolis parameter in the growth rate if all other parameters are 
unchanged, our experiments show no such preference. 


2 



3. Introduction 


Conditional instability of the second kind (CISK) (Ooyama 1964 and Chamey and 
Eliassen 1964; hereafter CE), had till not long ago been a popular concept in interpreting 
tropical cyclogenesis. Its variant, the frictional wave-CISK, has been a contender in the 
recent years in interpreting the Madden-Julian oscillation (e.g., Wang 1988, Xie and 
Kubokawa 1990.) Recently these ideas have come under criticism. Foremost among the 
critics are Emanuel, Neelin, and Bretherton (1994). They pointed out that in CISK the 
external view of the convective heating is implied; that is convective heating is specified 
without the quasi-equilibrium consideration. Among the other critics, Craig and Gray 
(1996) demonstrated in a numerical model of tropical cyclone (TC) that CISK does not 
work as well as another contending theory, the wind-induced surface heat exchange 
(WISHE, Emanuel et al. 1994) theory, in the sense that the intensification rate was found 
to be relatively insensitive to the frictional drag coefficient. This paper does not intend to 
assess the recent developments related to CISK. Instead, it focuses on two other aspects 

of CISK. 

The CISK theory has two prominent features. One is the important iole assigned to 
surface friction. In the CISK theory the low-level convergence is induced by surface 
friction, implying that surface friction is necessary for TC. Also according to the CISK 
theory (Chamey and Eliassen 1964), surface friction plays a dual role in the sense that 
besides its damping role it, through the Fictionally induced low-level convergence, also 
brings in moisture and thus plays an energy enhancing role. The other prominent feature 
is the idea of Ekman pumping (which is another name for the frictionally-induced low- 
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level convergence) efficiency (Charney 1971, 1973), which yields a TC growth tate 
proportional to the Coriolis force. In this paper these two features will be examined by 
focusing on the sensitivities of TC to surface frictional coefficient and the Coriolis 
parameter using an axisymmetric cloud-resolving model. 

2. Experimental setup 

The model used is an axisymmetric version of the latest version of the Goddard 
cloud ensemble model (Tao et al. 2002). An earlier version of this model was 
documented in Tao and Simpson (1993). The model variables include 3d velocity, 
potential temperature, perturbation pressure, turbulent kinetic energy and mixing ratios of 
all three water phases. An important feature of this type of models is that, since 
individual clouds are resolved, no cumulus parameterization is used. Therefore, unlike 
the meso-scale models, the success of this model does not depend on cumulus 
parameterization. However, the realism of the model still depends on the quality of the 
microphysics and turbulence parameterizations. Also the limitations of the 2-d structure, 
such as the lack of interaction between (azimuthal) mean flow and eddies, should be 
acknowledged. The model has 33 stretched vertical grids, which vary from 30 m near the 
bottom to 1140 m near the top. The vertical domain is about 21 km. The horizontal 
domain covers a radius of 924 km (770 grids) with a grid size of 1.25 km. The time step 
is 10 seconds. An open lateral boundary condition is used. The sea surface temperature 
is specified at 29°C. Unless specified otherwise, the Coriolis parameter is set at its value 

at 15°N. The surface fluxes are computed by the bulk scheme developed by Liu et al. 
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(1979). The Newtonian cooling is used with a time scale of 12 hours. The reference 
vertical temperature profile used in the Newtonian cooling is that of Jordan (1958). The 
Newtonian cooling models not only the effect of radiative processes (which would 
require the use of radiative equilibrium temperature profile as the reference profile 
toward which the temperature profile is relaxed to) but also the effect of processes with 
scale larger than that of TC. Our experiments showed that if the Newtonian cooling is 
replaced by the radiation package that exists in the original 2D model (short wave 
parameterization developed by Chou and Suarez (1999) and long wave parameterization 
developed by Chou et al. (1995), Chou and Kouvaris (1991), Chou et al. (1999) and 
Kratz et al. (1998)) that simulates radiative processes only (see Tao et al. 2002 for 
description, the simulated TC is slightly stronger. Other authors have also obtained 
reasonable simulation of TC with full radiation package (e. g„ Craig 1996). However, 
the use of full radiation package doubles the use of computer time for the model, we thus 
decided to proceed with the TC simulation with the Newtonian cooling. 

The initial condition has a temperature uniform in the radial direction and its 
vertical profile is that of Jordan (1958). The initial winds are a weak vortex in the 
tangential winds superimposed on a resting atmosphere. That is. 

r ^ 

v = max(0., v ), v, = (^) r exp(l. )(!•-—) 

r r 

max) max 1 

where v max = 10 m/s, r max =105 km, and z T =16 km. The initial sea level pressure is set at 
1002 hPs everywhere. Since there is no surface pressure perturbation, the initial 
condition is not in gradient balance. This initial cond.tion (along with the other model 
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settings) is enough to lead to cyclogenesis, but the initial stage (say, the first 24 his) is not 
considered to resemble anything observed. An initial condition that is in gradient wind 
balance would have a surge of upward motion at the center due to the Ekman pumping at 
the start of the integration (as surface friction reduces the Coriolis force, the difference 
between the pressure gradient force and the sum of the Coriolis force and the centrefugal 
force pushes the air in the boundary layer toward the center), which does not resemble 
anything observed either. A more satisfactory initial condition should have a balance 
among pressure gradient force, the Coriolis force, the centrifugal force, and surface 
friction; but that is more difficult to specify. Our initial condition avoids the initial 
artificial upward motion associated with the balanced gradient wind initial condition; thus 
the possibility of genesis being set off by the initial condition is avoided. 

The control experiment reaches quasi-equilibrium after about 10 days with an 
eye-wall extending from 15 to 35 km radius. Fig. 1 shows the results averaged over the 
24-hour period ending day 15 hr 0, which are reasonable in comparison with the 
observations. The maximum tangential wind of more than 65 m/s is located at about 20 
km radius at the top of the boundary layer. The minimum sea level pressure is about 955 
hPa. The maximum precipitation rate reaches more than 1 10 mm/hr. The warm core 
temperature perturbation reaches more than 12°C. Thus the overall simulation is deemed 

successful and the model is considered adequate for our experiments. None of our 
experiments has shown eye-wall replacement. 

3. Sensitivity experiments 
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Sensitivity to surface friction 


Three experiments identical to the control but with zero, half and double surface 
drag coefficient were done. The experiment without surface friction has also obtained the 
cyclogenesis process. Thus, surface friction is not a necessary condition for cyclogenesis 
or for the mature TC. However, without surface friction the TC has a very unrealistic 
structure. Fig. 2 shows the results averaged over the 24-hour period ending at day 15 hr 
0. The eye-wall in this case is quite far (about 70 km) from the center. The maximum 
tangential wind reaches 50 m/s at around 70km from the center. This maximum wind 
speed is weaker than that in the control experiment. However, since it occurs at a larger 
radius, it provides more kinetic energy. Unlike the control experiment, the maximum 
tangential wind is not very close to the surface. The converging flow toward the eye-wall 
obviously can exist without the presence of surface friction. However this converging 
flow is not confined in the low-levels. This experiment supports one of the central ideas 
of CISK that the low-level convergence in the tropical cyclone is induced by surface 
friction. In other words, surface friction is necessary for the existence of low-level 
convergence. Moreover, surface friction is important in giving the TC its observed 
smaller size and weaker intensity. The sensitivity of the low-level convergence to surface 
friction will be examined shortly. Without surface friction to reduce surface wind speed 
the converging air cannot reach very close to the center, resulting in an eye-wall with a 
large radius. This is due to the fact that the converging low-level air, in the absence of 
surface friction, gains high enough tangential wind at large radius to yield high outward 
centrifugal force and Coriolis force to prevent air from getting close to the center (l. e„ to 
counter the pressure gradient force.) 
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In the experiments with double and half surface friction coefficients the results 
show that higher surface friction renders a smaller and weaker TC bu, wt.h earlier and 
more intense cyclogenes, s process in terms of minimal sea level pressure (Fig. 3). Since 
the sea level pressure drop at the TC center depends on the cancellation between the 
upper-level divergence and the low-level convergence and, as will be soon shown, higher 
surface friction leads to higher low-level convergence and the more intense cyclogenes, s 
process helps evacuate air from the TC center leading to a more rapid sea level pressure 
drop. Similar to what was obtained by Craig and Gray (1996) using a mod.fied Rotunno 
and Emanuel ( 1987) model, the growth rate of the trap, cal cyclogenesis processes is not 
very sens, live to the surface friction if the minimum sea level pressure is used as a 
measure. Fig. 3a shows that the minimum sea level pressure drops at about the same rate 
for ail the cases. However, if the growth rate is measured in terms of total kinetic energy 
in the model domain, there is sensitivity to the surface friction. Fig. 3b shows that the 
growth rate of total kinetic energy, divided by the horizontal area of the model, is higher 
when surface friction is smaller. Specifically, all three (control, half and double) 
experiments reach about the same amount kinetic energy at hr 120 and the onset time is 
latest for the lowest surface friction case. 

Thus far we have shown that surface friction is important for low-level 
convergence and for a realistic simulation of TC. To further study the role of surface 
friction in TC, we need more experiments. In another senes of expenments the results of 
the control run at hrs 48, 96 and 240 (the incipient, growing, and mature phases) are used 
as the initial conditions for integralions of one day duration with half and double surface 
drag coefficient. Fig. 4.a shows the minimum sea level pressure of these expenments. It 
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reveals that in the 24 hour period starting from hr 48 (the incipient stage) doubling 
surface friction leads to a more rapid decrease of the minimum sea level pressure. Fig. 
5 .a shows the 24 hr averaged boundary layer (bottom 1 km) mass convergence for the 
control, half and double surface friction cases as functions of radius for these 
experiments. Apparently, higher surface drag increases the low-level convergence in all 
phases of the TC development. The increase in the low-level convergence is 
accompanied by an increase in the precipitation rate (Fig. S.b). In the 24 hour period 
starting hr 48 increasing surface friction leads to a faster minimum sea level pressure 
drop and a faster growth of total kinetic energy (Fig. 4.b). In the 24 hour periods starting 
from hr 96 (in the rapid intensification period) and hr 240 (in the mature stage) increasing 
surface drag also leads to an increase in low-level convergence but a slower drop of the 
central sea level pressure. The low-level convergence increase is the result of the 
reduction of the azimuthal wind and thus the reduction of the radially outward Coriolis 
force and centrifugal force to oppose the pressure gradient force. In the same periods 


increasing surface drag has decreased the total kinetic energy. 

These experiments reveal that the response of the cyclone to surface drag 
changes, in terms of minimum sea level pressure and total kinetic energy, \anes 
according to the developmental stage of the cyclone; but the response in the low level 
convergence does not vary. Thus the Ekman pumping concept, which states that low- 
level convergence is related to surface friction and which implies that low-level 
convergence increases when the drag coefficient is increased, appears to be valid. 
However, increased Ekman pumping yields faster TC growth only in the incipient phase. 
This is consistent with the finding that higher surface friction leads to earlier onset of the 
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genesis process (Fig. 3.a). Although higher surface friction can lead to earlier onset in 
terms of sea level pressure, it leads to higher area averaged kinetic energy only in the 
short early hours (Fig. 3.b). Therefore surface friction does play a dual role (the 
traditional energy damping role and the role of inducing moisture convergence in the 
boundary layer which supply latent energy to the core). However, with the minor 
exception in the early hours the net effect of higher surface friction is still to remove 
more kinetic energy from the system. 


Sensitivity to the Coriolis parameter 

Fig. 6 shows the minimum sea level pressure as a function of time for 
experiments with various values of the Coriolis parameter corresponding to various 
latitudes (while keeping everything else the same.) It shows that all expenments reach a 
minimal sea level pressure of about 930 hPa, but the onset is earlier for lower Coriolis 
parameter. The peak growth rate tends to become smaller with increasing Coriolis 
parameter. However, the minimum sea level pressure is not necessarily the only measure 
of the TC intensity. Fig. 7 shows the total kinetic energy, divided by the area of the 
model horizontal domain, of these experiments. It also shows that higher Coriolis 
parameter gives later onset. The maximum total kinetic energy that is attained depends 
on the Coriolis parameter. For values of the Coriolis parameter that corresponds to 
latitudes lower than 30N, the higher the Coriolis parameter the larger the maximum 
kinetic energy. For other values of the Coriolis parameter the maximum total kinetic 
energy does not very much depend on the Coriolis parameter. The maximum growth rate 
in terms of total kinetic energy increase appears to occur for the case of Coriolis 
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parameter corresponding to that of 35N. This figure also shows that the rates of total 
kinetic energy growth during the period of fastest growth are not very different between 
35° and 60° of latitude. At 90° the TC takes longer to reach its maximum total kinetic 

energy than at the latitudes between 35 and 60. Fig. 8 shows the precipitation rate 
averaged over day 10 to day 15 for various f values (their corresponding latitudes). The 
most optimal latitude for precipitation appears to be 20N in terms of total precipitation, 
eye-wall size and intensity. Below 30N the radius of the eye-wall increases as the 
latitude is raised (in computing 0: whereas above 30N there is no strong dependency. 
These results do not support the Ekman pumping efficiency idea of Chantey (1973) m 
interpreting the TC growth, which yields an optimal latitude at the poles in terms of 
growth rate for these experiments. The Coriolis parameter has two effects on the tropical 
cyclone. First, it provides inertial stability, which suppress convection and therefore the 
TC itself. This can be understood through the equivalence between f and vertical 
stability (Chao and Chen 2001). Second, the Coriolis parameter also generates azimuthal 
wind, which enhances the surface wind speed and thus evaporation. These two effects 
oppose each other in determining the optimal latitude. The former attracts the TC 
towards the equator and the latter towards the higher latitudes (i.e. the poles). 

4. Summary 

In summary, in this work we have examined the sens, unities of trop.cal cyclones 
in an axisytnmetic cloud-resolving model to surface friction and the Coriolis parameter. 
Although the tropical cyclogenesis process can occur without surface 
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developed TC is very unrealistic. Without surface friction the radius of the eye-wall is as 
large as 70 km and the maximum precipitation rate is more than one order of magnitude 
less. The total kinetic energy is much greater but the minimum surface pressure is much 
higher. Also the strongest inflow is not at the low-levels. Surface friction gives the TC 
the observed smaller size, and lower intensity. It also confines the inflow in the boundary 
layer. Further experiments with varying surface friction coefficients reveal that with 
higher surface friction coefficient cyclogenesis occurs sooner and faster but it terminates 
sooner also, giving a weaker TC. In all phases of the cyclogenesis process higher surface 
friction giving higher boundary layer convergence. However the responses in minimum 
sea level pressure and total kinetic energy are different in different phases. In the very 
beginning phase of the cyclogenesis, higher surface friction leads to faster drop ol the 
minimum sea-level pressure and an increase in total kinetic energy, just the opposite to 
what happens in the later phases. Our study also shows that contrary to the Ekman 
pumping efficiency idea in the interpretation of TC growth, which prefers the highest 
value for f, tropical cyclogenesis process has a different preferred value for f, its value at 
35N if measured in terms of fastest rate of drop of the minimum sea level pressure and in 
terms of fastest growth of total kinetic energy that is attained during the cyclogenesis 

processes (Figs. 6 and 7). 
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Figure Caption 


Figure 1 . Results of the control experiment averaged over the 24 hr period ending at hour 
0 of day 15. (a) tangential wind (m/s), (b) radial wind (m/s), (c) vertical 
velocity (10" 2 m/s), (d) temperature minus reference temperature (°C), 

which a function of pressure only, (e) angular momentum (m2s-l), (f) 
potential vorticity (10~^ m 2 kg~^ s*^) as functions of height and radius, (g) 
sea level pressure (hPa) as a function of radius, and (h) precipitation 
(mm/day) as a function of radius. 

Figure 2. Same as Figure 1 but for the experiment without surface friction. 

Fieure 3. (a) Minimum sea level pressure (in hPa) as functions of time for control, half, 
double and no surface friction. 

(b) same as (a) but for total kinetic energy per unit area (in 10 4 J/m-.) 

Figure 4. (a) Minimum sea level pressure (in hPa) as functions of time for the control 
experiment and for the double and half surface friction experiments for 
three 24 hour periods starting from the control and (b) same as (a) but for 
total kinetic energy per unit area (in 10 4 J/m 2 .) 

Figure 5. a. Mass flux in radial direction in the bottom one kilometer in 10 9 kg/s as 
functions of radial distance from the TC center for the control and half and 
double surface friction averaged over the 24 hour period starting (a) houi 
48, (b) hour 96 and (c) hour 240. 

Figure 5.b. Same as Fig. 5. a but for precipitation rate (in mm/hr). 
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Figure 6. 

Figure 7. 
Figure 8 


Minimum sea level pressure (in hPa) as a function of time for various values of 
the Coriolis parameter, labeled by the corresponding latitude 
Same as Fig. 6.a but for total kinetic energy per unit area (in 10 4 i/m 2 .) 

Same as Fig. 6 but for precipitation rate (mm/day) averaged between day 10 

and day 15. 
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